

library(forecast)
library(tseries)   #调入两个库

data <- read.csv("D:/yanyi/data/qq.csv")  #读取csv文件

ax=ts(data$co2,start=1820)   #以年份和co2浓度构成时间序列

#画出原始数据图像
plot(ax,col=4,pch=8,type="o",lwd=2)


#使用Holt双参数指数平滑模型，α和β的选择根据时间序列自适应
x.fit=HoltWinters(ax,gamma=F)
x.fit

#画出样本内预测图像
plot(x.fit)


#对未来30年进行预测
x.fore=forecast(x.fit,h=30)
#输出预测结果
x.fore

#画出未来30年预测图像
plot(x.fore,col=2,lwd=2)
